function [ fout, tout ] = get_blip_python( amiga_model, ... 
    audio_filter_mode, unquatized, kind, mode )

fs = get_sampling_fq();

if (nargin < 3)
    unquatized = true;
end

if (nargin < 4)
    kind = 'lsim';
end

if (nargin < 5)
    mode = 'pal';
end


if strcmpi(kind, 'lsim')
    if unquatized
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-unquantized-lsim-%s-%s.csv', amiga_model, mode));
    else
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-lsim-%s-%s.csv', amiga_model, mode));
    end
elseif strcmpi(kind, 'org')
    if unquatized
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-unquantized-org-%s-%s.csv', amiga_model, mode));
    else
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-org-%s-%s.csv', amiga_model, mode));
    end
elseif strcmpi(kind, 'bq')
    if unquatized
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-unquantized-bq-%s-%s.csv', amiga_model, mode));
    else
        csv = importdata(...
            sprintf('..\\..\\genblep\\bleps-bq-%s-%s.csv', amiga_model, mode));
    end
else
    error('bad lsim param');
end

data = csv.data;

switch audio_filter_mode
    case 'led filter off'    
        index = 2;
    case 'led filter on'
        index = 3;
    case 'unfiltered'
        index = 4;
    otherwise
        error('bad audio_filter_mode parameter');
end

fout = data(:,index)';
tout = linspace(0, (size(fout, 2)-1) / fs, size(fout, 2));

% step 1->0 switch to 0->1
fout = fout(1) - fout;
fout = fout / fout(end);

% blep -> blip
fout = destep_calc(fout);

end